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Production and decay of ft® baryons is studied with ~ 230 fb -1 of data recorded with the BABAR 
detector at the PEP-II e + e~ asymmetric-energy storage ring at SLAC. The ft® is reconstructed 
through its decays into i?~7r + , i?~7r + 7r - 7r + , and S~K~tt + 7t + final states. 

The invariant mass spectra are presented and the signal yields are extracted. Ratios of branching 
c~| ! fractions are measured relative to the ft® — > i?~7r + mode 

B(ft° — > K~it + it + ) 

B{ft{-^ft-^) = °- 31 ± ai5(stat ° ± °- 04 ^ st -)' 
S 1 B(ft° ^ ft-ir + ir-ir+) , M „ . 

bU~o-*+) < °- m (90%CL) - 

The momentum spectrum (not corrected for efficiency) of ft® baryons is extracted from decays 
into i7~7r + , establishing the first observation of ft® production from B decays. 
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1 INTRODUCTION 



The f2® (ess) is a (J p = \ )^ ground state baryon with a mass of m^o = (2697.5 ± 2.6) MeV/c 2 
and a lifetime of Tqo = (69 ± 12) fs [1]. Since the first evidence for f2® production and decay in 
1984 [2] in the decay mode Q° c -► E~K~ ~tt + tt + , the f2® baryon has been seen by a number of 
different experiments [3-8] in various decay modes, with strong evidence reported in the decays 
J2° -> ^"vr + , -> ~ , "K"vr + 7r+, and 12° -> S + K~K-tt + . To date, only one 5cr observation 
has been reported by a single experiment, combining two modes [8]. No observation of a single 
exclusive f2® decay mode at the 5<r level has been reported and only a few of its decay modes have 
been observed. The measurements of the ratios of branching fractions still have large uncertainties 
and the production mechanisms of the remain largely unexplored. 

The large amount of data collected at the B factories allows a more detailed analysis of the 
properties of the . In this paper, a study of the through the decay modes 4 l?~7r + , J2~7r + 7r - 7r + , 
and £~ K~ir + ir + is described. 

The invariant mass spectra are presented and the ratios of branching fractions relative to the 
J? - 7T + decay mode are calculated. To study the production mechanism, the f2® momentum spec- 
trum (not corrected for efficiency) in the e + e~ rest frame is presented. 

2 THE BABAR DETECTOR AND DATASET 

The BABAR detector operating at the PEP-II e + e~ asymmetric-energy storage ring at SLAC con- 
sists of a tracking system for the detection of charged particles, a detector of internally reflected 
Cherenkov light (DIRC), an electromagnetic calorimeter (EMC), and an instrumented flux return 
(IFR). The tracking system, contained in a 1.5-T magnetic field provided by a superconducting 
solenoidal coil, includes a 5-layer, double-sided silicon vertex tracker (SVT) and a 40-layer drift 
chamber (DCH). The EMC consists of 6580 CsI(Tl) crystals. Information from the DIRC and 
the energy loss information from the SVT and the DCH are used for charged particle identification 
(PID). The IFR is segmented and instrumented with resistive plate chambers. The BABAR detector 
is described in detail elsewhere [9]. 

This analysis is based on data taken at the T(4S) resonance and ~ 40MeV/c 2 below, referred to 
as on-peak and off-peak data, respectively. The integrated luminosity of the data used corresponds 
to ~ 225 fb -1 in each of the decay modes involving an fi~ , and ~ 230 fb" 1 in the 5~ decay 
mode. 

To minimize the selection bias, all selection criteria are optimized on Monte Carlo simulated 
event samples at least as large as the data sample. Samples of e + e~ — > qq (q € {u, d, s, c}) events 
are used to study the background. Monte Carlo samples are used to study the signal properties 
and to evaluate the selection efficiencies. The signal samples are generated with JETSET [10], 
assuming uniform phase space for the 4?° decays. 

3 EVENT SELECTION 

The fi~ ("~) is reconstructed in its AK~ (^tt - ) final state. Candidates for A decays are re- 
constructed in their pir~ final state. The branching fractions [1] of the intermediate hyperons are 
B{Q~ AK~) = (67.8 ± 0.7)%, B{S~ -► Air~) ~ 100% , and B(A -► pvr~) = (63.9 ± 0.5)%. 

*The quantum numbers have not been measured, but are assigned in accord with the quark model. 
4 Sinmltaneous treatment of the charge conjugate mode is always implied throughout the note. 
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3.1 The Hyperon Selection 

All hyperons in this analysis (A, S~ , and are long-lived particles with a typical decay length 
of several cm in BABAR. Each hyperon is identified by reconstructing its decay vertex, which is 
required to be clearly displaced from that of the parent particle. In the selection of each of the 
intermediate hyperons, only candidates with an invariant mass within 3a of the central value are 
selected, where a denotes the invariant mass resolution. These candidates are then subjected to a 
kinematic fit, constraining the mass of the candidate to its nominal value. 

Candidates of A baryons are formed from a pair of tracks of opposite charge, where the positively 
charged track must satisfy the PID requirements for a proton. Each A candidate is combined with 
a 7r~ (K~) to form a E~ (&~) candidate. The K~ candidate track from the Q~ decay must 
satisfy the PID requirements for kaons. 

Candidates for fi® decays for each final state are formed by combining the reconstructed fi~ 
and S~ baryons with the required number of mesons (7r _ , 7r + , or K~). Throughout this document, 
mesons that are the direct daughters of the candidate are referred to as 'primary mesons'. 

3.2 The i?° Candidate Selection 

The kinematics of the signal decays and the backgrounds are different for each of the i?° decay 
channels considered in this analysis. The Q-values, which are a measure of the available kinetic 
energy of the decay products, are 882 MeV, 603 MeV, and 595 MeV for the I? - 7r + , i?~7r + 7r~7r + , and 
S~ K~ir + ir + decay modes, respectively. Therefore, the selection criteria are optimized for each 
mode separately. Particle identification is required for each of the primary mesons. The selection 
criteria described below are then applied to the candidates in order to suppress backgrounds. In 
addition, a common minimum p* of 2.8GeV/c is required for candidates used in the measurement 
of the ratios of branching fractions, where p* is the £2® momentum in the e + e~ rest frame. This 
is above the kinematic limit for production from B decays {p* = 2.02GeV/c). 

• -> [2-TT+: 

The transverse flight length of the ft~ candidate in the xy-plane calculated with respect to the 
event vertex is required to exceed 2 mm. The signed 5 flight length of the A, measured from 
the fi~ decay point, must exceed 1.5 mm. For the primary pion, a minimum momentum 
of 200MeV/c in the laboratory frame is required. A minimum of 0.1% is required for the 
X 2 -probability of the kinematic fit for each intermediate hyperon. 

• Q° c -> Q-TT+TT-TT+: 

The transverse flight length of the Q~ candidate in the xy-plane calculated with respect to 
the event vertex is required to exceed 2.5 mm. The signed flight length of the A, measured 
from the Q~ decay point, must exceed 2 mm. The vector sum of the momenta of the 7r + 7r~7r + 
system in the lab-frame is required to exceed 650MeV/c. A minimum of 0.1% is required for 
the x 2 -P r °bability of the kinematic fit for each intermediate hyperon. 

• Q° c -> ~-K-tt + tt + : 

The minimum requirement on the flight length of the S~ is 4.5 mm and the signed flight 
length of the A with respect to the S~ decay vertex is required to be larger than zero. A 
minimum x 2 probability of 10 -4 is required for the kinematic fit of the full decay chain. 

6 The signed flight length is defined as the dot product of the displacement vector (va — fn) of the A and the 
momentum vector of the A, where r denotes the 3D position of the vertex. 
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4 THE INVARIANT MASS SPECTRA 



The signal yield for each decay mode is extracted from an unbinned maximum likelihood fit to the 
invariant mass spectrum of the f2® candidates. In the fit, the signal lineshape is described by a sum 
of three Gaussians with a common mean, and the background is described by a first order 
mode) or second order (all other modes) polynomial. The widths and relative contributions of each 
Gaussian are fixed to the values obtained from samples of simulated f2® baryons produced in cc 
continuum and decayed in the considered mode. The mean mass, the yield of the signal peak, and 
the parameters in the polynomial for the description of the background are left free in the fit to the 
data. In the fit to the .f2~7r + 7r~7r + mode, due to the small signal, the mean of the signal lineshape 
is fixed to the weighted mean 6 (fi = 2694.6 MeV/c 2 ) obtained from the fits to the other two decay 
modes. 

The invariant mass spectra are displayed in Figures l(a-c) for the ,!?~7r + , Q~ tt + ir~ ir + , and 
E~K~ir + ir + modes. The results of the fits to the data, the selection efficiencies 7 , and the \ 2 



6 This is not necessarily the true mass of the J?° baryon. No studies of systematic biases due to energy loss 
corrections of the decay products or due to the uncertainties on the energy and momentum scale are carried out. 



7 A11 selection efficiencies quoted are for decays with p* > 2.8 GeV/c 2 
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Figure 1: Invariant mass spectra for fi® decays into (a) J7~7r + , (b) i?~7r + 7r~7r + , and (c) S~ K~tt + it + final 
states with p* > 2.8 GeV/c. The dots are the data. The binning is chosen according to the Half- Width-Half- 
Maximum of the signal peak obtained from signal Monte Carlo. The result of the fit is overlaid 
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Table 1 : Results from the fits to the invariant mass spectra in data. The yields (not corrected for efficiency) 
are given for each mode individually, as well as the selection efficiencies. In addition the x 2 probability for 
the fit, calculated in the mass window 2.6 < m < 2.8GeV/c 2 , and the significance S of the signal is quoted 
for the j?~7r + 7r~7r + and the S~ K~tt + tt + modes. 



Decay Mode 


Signal Yield 


Efficiency (%) 


prob(x 2 ) 


S 




138.5 ± 14.8 


8.35 ±0.07 


0.54 


17.8 




11.8 ±7.5 


4.41 ±0.09 


0.73 


2.4 




29.9 ± 13.6 


5.63 ±0.10 


0.85 


3.4 



probability for each fit are summarized in Table 1. In addition, the significance S for the signal, 
denned as y/2A log C, calculated from the difference in the log-likelihood (log C) for a fit with and 
without a signal lineshape, is included. 

5 PRODUCTION MECHANISM FOR Q Q C BARYONS 

Insight into the production mechanism for baryons is obtained from the p* spectrum. For this 
study, the fi® — ► I? - 7r + decay mode, which has the largest signal yield of all modes in this analysis, 
is used. The signal yield as a function of p* is measured up to 4.4GeV/c in eleven intervals, each 
400MeV/c wide. For each interval, the f2~ir + invariant mass spectrum is fit with the lineshape for 
the signal fixed to that obtained from the full signal Monte Carlo sample. A polynomial is used to 
describe the background. No significant variation in the central value of the signal peak is observed 
as a function of p* . It is therefore fixed to the common mean for all intervals in the fit. 

The measured signal yields obtained from the combined on-peak and off peak data sets can be 
compared with those from the off-peak data set, displayed in Figures 2(a) and 2(b), respectively. 
The dots represent the data and the error bars correspond to the statistical uncertainty only. The 
solid horizontal bars correspond to the predicted spectrum for production from cc continuum 
Monte Carlo; the thickness of the bars correspond to the statistical uncertainty. No correction for 
selection efficiency is applied to either distribution. The corresponding distribution from signal 
Monte Carlo is normalized such that its integral corresponds to that in data for p* > 2.5GeV/c in 
Figure 2(a), and for the full p* range in Figure 2(b). 

A clear two-peak structure is evident in the distribution from the combined on-peak and off- 
peak data sets. The peak at high p* is consistent with f2® production as predicted from continuum 
signal Monte Carlo and the off-peak data. The p* spectra from data and Monte Carlo show good 
agreement within the experimental uncertainties in this region. The peak in the p* region below 
2.02GeV/c provides clear first evidence for production from B decays. This interpretation is 
substantiated by the absence of the corresponding peak in the spectrum extracted from off-peak 
data only, taken below the BB threshold. 

6 PHYSICS RESULTS 

6.1 Ratios of Branching Fractions 

The yields for ft® signal events, extracted from the invariant mass spectra in the data and corrected 
for selection efficiency and acceptance effects. These are extracted from reconstructed signal Monte 
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Figure 2: The signal yield as a function of the p* of candidates (a) from the combined on-peak and off- 
peak data sets and (b) from off-peak data only. The dots are the data and the vertical error bars correspond 
to the statistical uncertainty only. The solid horizontal bars correspond to the predicted distribution for J7° 
production from cc continuum signal Monte Carlo. The thickness of the bars correspond to the statistical 
uncertainty of the Monte Carlo sample. No correction for selection efficiency is applied to any of the 
distributions shown. The distributions are normalized to the same area for (a) p* > 2.5 GeV/c and (b) the 
full range. In (a), clear evidence for f2® production from B decays is visible at lowp*. ft is absent in off-peak 
data (b), collected below the BB production threshold. 



Carlo events that pass the same selection criteria as events in data. The efficiency corrected yields 
are then used to calculate the ratios of branching fractions relative to the — * f2~ir + mode, 
yielding 

B(f2° c -> ^-^-vr+) = lg ± 1Q(gtat ± a03(sygt)) 



B{Q Q C -► J7-7T+) 
B{Q° C -» ~-K-it + tt + ) 



0.31 ± 0.15(stat.) ± 0.04(sys.). 



B{Q° C -> f2-Tr+) 

The individual contributions to the systematic uncertainties are discussed in detail in Section 6.2. 



6.2 Systematic Studies 

Although the decay topologies are slightly different in the various decay modes, the systematic 
uncertainties on the hyperon selection efficiencies largely cancel in the ratios of branching fractions. 
The following sources of systematic uncertainties are considered in the measurement of the ratios 
of branching fractions and summarized in Table 2. 

• Monte Carlo Simulation: 

The statistical uncertainty on the selection efficiencies is taken into account in the systematic 
uncertainty. 

In multi-body decay modes, the decay can occur via intermediate resonances. Therefore, sig- 
nal Monte Carlo samples containing the ='*(1530) A'"7r + , the S~K* q tt + , and the S*{1530) K*< 
decay modes are generated. The difference in selection efficiencies relative to the uniform- 
phase-space sample is included as an uncertainty on the selection efficiency. For the i?~7r + 7r~7r + 
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mode, the systematic uncertainty associated with intermediate resonances is assumed to be 
the same as that for the topologically similar mode S~ K~it + tt + . 

The p* spectrum of the f2® in data and Monte Carlo is slightly (but not significantly) different. 
The effect on the selection efficiency is calculated and a systematic uncertainty accounting for 
the difference is assigned. This contribution is correlated for the decay mode in the numer- 
ator and in the denominator and is treated accordingly in the calculation of the systematic 
uncertainty on the ratio of branching fractions. 

• Extraction of the Signal Yield: 

In the fit to the invariant mass spectra, the description of the background shape is varied for 
the purpose of estimating the systematic uncertainty. In addition, the width of the widest 
Gaussian in the description of the signal shape is varied by a factor of two, and the fit range 
is varied from 2.6 < m < 2.8GeV/c 2 to 2.4 < m < 3.0GeV/c 2 . The observed variations in the 
signal yield are added as a systematic uncertainty. 

• Particle Identification and Tracking: 

The number of primary mesons is different for some of the modes. A systematic uncertainty of 
1% is added for each additional primary pion, estimated from the uncertainty on the efficiency 
of the pion identification. A systematic uncertainty of 1% is assigned to account for different 
kaon identification algorithms used for the primary and non-primary kaons. 

In order to account for the difference in tracking efficiency in data and Monte Carlo, a 
correction of 0.25% with a systematic uncertainty of 1.4% is applied per track. In the ratio 
of branching fractions this amounts to a 0.5% correction to the selection efficiency with a 
systematic uncertainty of 2.8%. 

• Branching Fraction: 

The Q~ is identified in its decay to AK~, which has a branching fraction of (67.8 ± 0.7)%. 
This uncertainty in the branching fraction is added to the systematic uncertainty for the ratio 
involving the S~ only. The uncertainty in the A branching fraction to pn~ cancels in the 
ratio of branching fractions for all modes. 



Table 2: Systematic uncertainties considered in the measurement of the ratios of branching fractions. The 
individual contributions are given, and added in quadrature to determine the total systematic uncertainty. 
Dashes indicate sources that are assumed to cancel in the ratio of branching fractions. 







B(~-K-ir + ir + ) 
B(f2-ir+) 


Monte Carlo Statistics 


0.004 


0.006 


p* Rewcighting 


0.001 


0.004 


Resonance Structure 


0.003 


0.005 


Extraction of Signal Yield 


0.028 


0.032 


Particle ID & Tracking 


0.006 


0.011 


f2~ Branching Fraction 




0.003 


Multiple Candidates 


0.002 


0.004 


Total Systematic Uncertainty 


0.03 


0.04 
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• Multiple Candidates: 

A possible source of background comes from the presence of multiple fi® candidates in an 
event, which share one or more tracks. Predominantly, the same hyperon combines with one 
or more primary tracks to form such multiple candidates. In general, these candidates are 
distributed over a large mass range. However, in cases where these candidates share tracks, 
their masses might be correlated, which could lead to a pile-up of candidates in the signal 
region. The mass distribution of in events with multiple candidates is studied in data 
as well as in Monte Carlo signal samples. We select events in which one of the candidates 
falls inside a ±3 HWHM mass window around the nominal peak position. In data, two, four, 
and four events with such multiple candidates are observed in the Q~ir + , ,!?~7r + 7r~7r + , and 
E~ K~it + it + , respectively, where the second candidate also lies in the signal region defined 
above. From studies of these candidates in a larger mass window, these multiple candidates 
are observed to be evenly distributed in mass for all modes. Therefore, they form part of the 
background when the mass spectra are fit with a polynomial. This indicates that there is no 
statistically significant peaking under the signal peak in data. 

The relative sizes of the peaks from incorrectly reconstructed multiple candidates and true 
signal candidates in Monte Carlo are 0.2%, 1.2% and 1.3% in the f]~ir + , f2~ir + Tr~ir + , and 
E~ K~ir + ir + mode, respectively. This fraction is assigned as a systematic uncertainty. 

• Other Sources of Peaking Background: 

Possible sources of peaking backgrounds are studied with continuum Monte Carlo samples. 
The number of peaking background events observed is consistent with zero and therefore no 
systematic uncertainty is assigned. 

The individual sources of systematic uncertainties are added in quadrature. The total uncer- 
tainty on each ratio is given in Table 2. 

6.3 Limit Calculation for the Q° c -> !? tt+tt tt+ Mode 

No significant excess of signal events over the background is observed in the I2~7r + 7r~7r + mode. 
Therefore, a limit at the 90% confidence level (CL) on the ratio of branching fractions is calculated. 
The limit is obtained from a Monte Carlo calculation for the individual modes, using the measured 
yields and the statistical and systematic uncertainties as inputs. All uncertainties are assumed 
to be Gaussian. Integrating the positive part of the distribution of ratios of branching fractions 
obtained from these Monte Carlo experiments, an upper limit at the 90% confidence level of 



< 0.30 (90%CL) 



is obtained. 

6.4 Cross-checks 

Cross-checks for the measurements of the ratios of branching fractions are performed to verify the 
stability of the result. 

• Charge asymmetry: 
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To obtain the main result, no distinction is made between decays of the 17° and the anti- 
baryon 17 C . However, the selection efficiencies might be different due to differences in the 
interaction in material of particles and anti-particles. Particle and anti-particle candidates 
are selected and studied separately. The observed difference in efficiency-corrected yields for 
mode and anti-mode is largest in the 17° decay into S~A'~7r + 7r + , with a difference of 1.6a. 
All ratios of branching fractions are found to be consistent within the statistical uncertainties. 

• Choice of the p* range: 

The main results for the ratios of branching fractions are obtained for a minimum p* of 
2.8GeV/c. The ratios of branching fractions, however, are expected to be independent of 
that choice. As a cross-check, these ratios are recalculated, requiring a p* in the range 
2.6 GeV/c < p* < 3.2 GeV/c and p* > 3.2 GeV/c. All results are found to be consistent within 
statistical uncertainties. 



7 SUMMARY 

Data recorded with the BABAR detector are analyzed to study 17° production and decays. The 17° 
is reconstructed through its decays into and l?~7r + 7r~7r + (using 225 fb^ 1 of data), and into 

Z-K-Tt+Tt + (using 230 ftT 1 ). 

Based on the momentum spectrum of the 17° in the e + e~ rest frame in the full data set and 
in off-peak data only, the first observation of 17° production from B decays is reported. 

In the i? - 7r + decay mode, the 17° is observed with a statistical significance of over 17a. This 
constitutes the first observation of the 17° baryon above the 5a level. 

From the observed signal yields, corrected for efficiency and acceptance, the ratios of branching 
fractions relative to 17° — > l7~7r+ are measured to be 



B{Q° C ->■ fl-vr+) 

g(!7 c ° -> S-K-TT+TT+) 



0.16 ± 0.10(stat.) ± 0.03(syst.), 
0.31 ± 0.15(stat.) ± 0.04(sys.), 



where the first uncertainty is statistical and the second is systematic. Due to the limited statistical 
significance for the decay mode 17° — > l7~7r + 7r~7r + , an upper limit on the ratio of branching 
fractions, 

B((2° - n-ir+) K °- 3 ° ( 90%CL )> 

at the 90% confidence level is set. All results mark a considerable improvement, both in statistical 
and systematic uncertainties, over the current world averages [1]. 
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